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Optimal disorder for segregation in annealed small worlds 
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Abstract. We study a model for microscopic segregation in a homogeneous system of particles moving on a 
one-dimensional lattice. Particles tend to separate from each other, and evolution ceases when at least one 
empty site is found between any two particles. Motion is a mixture of diffusion to nearest-neighbour sites 
and long-range jumps, known as annealed small- world propagation. The long-range jump probability plays 
the role of the small-world disorder. We show that there is an optimal value of this probability, for which 
the segregation process is fastest. Moreover, above a critical probability, the time needed to reach a fully 
segregated state diverges for asymptotically large systems. These special values of the long-range jump 
probability depend crucially on the particle density. Our system is a novel example of the rare dynamical 
processes with critical behaviour at a finite value of the small-world disorder. 

PACS. 89.75.Fb Structures and organization in complex systems - 02.50.Ey Stochastic processes 
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Segregation -the spatial separation of different ele- 
ments in a multi-component system- is a paradigmatic 
form of complex behaviour, present in a wide variety of 
natural phenomena. It can occur both at the microscopic 
level, driven by repulsive interactions between individual 
components, and at the macroscopic level, as in the case 
of phase separation due to collective processes. In many 
instances, it leads to the spontaneous emergence of spatio- 
temporal order, in the form of coherent structures with 
non-trivial evolution. 

In physical systems, spatial segregation is conspicuous 
in the vicinity of phase transitions, as the result of the 
amplification of microscopic fluctuations 1 . An effect of 
similar origin is observed in some elementary chemical re- 
actions limited by diffusion. It is well known, for example, 
that the binary annihilation of two diffusing chemicals A 
and B, A + B — > 0, leads to the formation of single-species 
spatial domains just after the very first stages of the pro- 
cess, starting from a uniform mixture of the two species 

El 

In the realm of biology, segregation at the molecular 
and cellular level plays a prominent role in morphogene- 
sis, i.e. in the differentiation of the various tissues, organs, 
and functional parts of living organisms |3j. At the level 
of ecosystems, environmental heterogeneities may lead to 
population segregation, with the appearance of specific 
localized niches, when the relative fitness of competing 
species varies in space 0]. In human populations, the com- 
bination of cultural diversity with xenophobic feelings is 
able to disrupt society, generating hatred and aggression 
where cooperative behaviour should otherwise emerge [S]- 
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Homogeneous ecological and social systems, with no 
substantial differences between their individuals, often ex- 
hibit a different form of segregation, not necessarily as- 
sociated with the formation of spatial structures. Single 
individuals -or small groups of closely related individuals, 
such as families or clans- may choose to separate from 
each other as much as possible, all over the available spa- 
tial domain. This spreading is related to a more efficient 
exploitation of the resources associated with land, and has 
certainly been crucial to the geographical dissemination of 
biological species over the ages 0]. 



In this paper, we explore a physically-inspired model 
for this kind of segregation. We pay particular attention 
to the dynamical process that leads, from a random distri- 
bution of the population, to a stationary state where indi- 
viduals have mutually separated such that no contiguous 
neighbours are found. Individuals are allowed to move dif- 
fusively and, occasionally, they can perform long jumps to 
randomly selected sites. This form of transport, inspired 
in the connection structure of small- world networks, is dis- 
cussed in the next section. In Sect. we present our first 
numerical results, and show that the time needed to reach 
full segregation depends non-monotonically on the proba- 
bility of long jumps. Remarkably, in fact, the time is min- 
imal for an intermediate value of that probability. Section 
13 is devoted to develop an analytical approximation of 
our problem, which makes it possible to explain some of 
the numerical results. Further simulations are presented 
in Sect. where we validate predictions of the analyti- 
cal approximation and complete the study of our system. 
Results are summarized and discussed in the last section. 
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Fig. 1. Possible events associated with the motion of a frus- 
trated particle. The moving particles is shown in gray. Left col- 
umn: long-range jump events. Right column: diffusion events. 



1 Segregation in annealed small worlds 

We consider an ensemble of N particles distributed on a 
one-dimensional array of L sites (L > N) with periodic 
boundary conditions. The density is p = N/L. An exclu- 
sion principle holds so that, at any time, each site may be 
occupied by at most one particle. Particles tend to seg- 
regate, and a particle is said to be frustrated if at least 
one of its two nearest-neighbour sites is occupied. Frus- 
trated particles are allowed to move, in order to reach a 
non-frustrated state. Motion of frustrated particles is gov- 
erned by annealed small-world transport, as described in 
the following. 

Small-world networks have been defined to capture 
two essential properties of real social structures, namely, 
the small average distance between any two individuals 
(small- world effect), and their high clustering They 
have been used, as toy models, to demonstrate the role 
of the underlying interaction pattern in dynamical pro- 
cesses of socio-economical and ecological inspiration 
EE]- These networks are built starting from an ordered 
lattice of relatively high connectivity. A prescribed frac- 
tion of links is reconnected at random, thus introducing 
quenched disorder and creating shortcuts between distant 
parts of the original lattice. Annealed small worlds consti- 
tute a variation of the same model, introduced for systems 
where propagation processes occur between sites [1011111 
Ej. In this version of the small-world model, propagation 
takes place between neighbour sites of the ordered lattice 
but, with a certain probability, jumps to randomly chosen 
sites may take place. Long-range jumps play the role of 
the shortcuts of the quenched original version. By anal- 
ogy, these random jumps are associated with disorder in 
the propagation process. 

In our model of segregation, a particle is chosen at ran- 
dom at each time step. If it is not frustrated, it remains 
at its site. If, on the other hand, the selected particle is 
frustrated, it is allowed to jump to a randomly selected 
empty site, anywhere in the lattice, with probability p. 
With the complementary probability, 1 — p, the motion 
is diffusive. If only one of the nearest-neighbour sites is 
occupied, the particle moves to the empty neighbour. If, 
otherwise, the two nearest-neighbour sites are occupied, 
the exclusion principle forbids motion. The possible in- 
stances are illustrated in Fig. ^ As discussed above, the 
parameter p can be associated with the level of disorder of 
the annealed small world. The duration of each time step, 



At = 1/N, is fixed so that, on the average, each particle 
is chosen once per time unit. 

We may assert a priori that the system has two well- 
differentiated regimes, depending on how the density p 
compares with 1/2. For p < 1/2, collective states where 
no particle is frustrated exist. For finite L, moreover, the 
population will necessarily reach one of such fully non- 
frustrated states, from any initial condition, in a finite 
time T. From then on, it will remain in a frozen configura- 
tion where all particles are mutually separated by at least 
one empty site. On the other hand, for p > 1/2, there will 
always be frustrated particles, as there is not enough place 
to accommodate the population in a fully non-frustrated 
state. As time elapses, the number of frustrated particles 
will eventually fluctuate around a well-defined stationary 
value. 

Following this discussion, our analysis of the present 
model focuses on the study of the evolution of the number 
of frustrated particles. In the low-density regime, where 
the population is able to accomplish full segregation, we 
record the average time T needed to reach the fully non- 
frustrated state, starting from a random distribution of 
particles. The average is performed over different real- 
izations of the initial state and of the stochastic evolu- 
tion. For high densities, on the other hand, we analyze 
the asymptotic number of frustrated particles. For conve- 
nience both in the numerical and in the analytical calcu- 
lations, instead of recording the number of frustrated par- 
ticles we follow the evolution of the number P of particle 
pairs, defined as the number of particles whose nearest- 
neighbour sites to the right are occupied. The time average 
of the pair density p = P/N acts as an order parameter for 
the transition between the two regimes. Our main param- 
eters are the particle density p and the jump probability 
p. 



2 Preliminary numerical results 

We begin by exploring the scaling properties of our model, 
for fixed p and p, as a function of the system size L - 
or, equivalently, of the number of particles N = pL. It 
results that, for sufficiently large systems, the evolution of 
the pair density p becomes size-independent when time is 
suitably scaled by the step duration At = 1/N. In further 
simulations, we fix L = 5 x 10 4 , for which this asymptotic 
independence of size already holds. 

Figures [3 and show the evolution of the pair density 
p for several values of the particle density p and the two 
extreme values of small- world disorder, p — and 1, re- 
spectively. Each curve has been obtained by averaging over 
100 realizations. For p — 0, as argued above, p vanishes at 
finite times if p < 1/2 and reaches a slightly fluctuating 
stationary level if p > 1/2. The two regimes are separated 
by the case p = 1/2, where p displays a power-law decay 
with an exponent close to —1/4. 

For p = 1, rather unexpectedly, the transition between 
the two regimes does not take place at p = 1/2 but at 
a lower density p c , somewhere between 0.3 and 0.4. Even 
when we know that, for p c < p < 1/2, fully non-frustrated 
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Fig. 2. Evolution of the pair density p for p = and several 
values of the particle density p. Curves without labels corre- 
spond, from bottom to top, to p — 0.6, 0.7, 0.8, and 0.9. The 
dotted line has slope —1/4. 
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Fig. 3. Evolution of the pair density p for p = 1 and several 
values of the particle density p. Curves without labels corre- 
spond, from bottom to top, to p = 0.6, 0.7, 0.8, and 0.9. 



states certainly exist, the system is not able to find any of 
them and the pair density approaches a stationary non- 
zero level. It must be recalled that, as discussed above, a 
finite population will always reach a fully non-frustrated 
state due to the effect of fluctuations. For the system size 
of our simulations, however, it is apparent from Fig.[3]that 
the associated time grows beyond our numerical reach as 
soon as p becomes larger than p c . 

To reveal this feature in more detail and, in particular, 
to analyze the effect of the small- world disorder, we study 
the time T needed to reach a fully non-frustrated state. 
Figure 01 shows the average of T, taken over 500 to 10 4 



Fig. 4. Average time T needed to reach a fully non-frustrated 
state as a function of the disorder p, for various values of the 
density p. Data without labels correspond, from bottom to top, 
to p = 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, and 0.325. For clarity, data 
have been joined with dotted lines. 



numerical realizations, as a function of the disorder p and 
for several values of the density p < 1/2. In agreement 
with our previous results, we find that, for p > 0.35, there 
is a critical level of disorder p c at which the average time T 
becomes very large. For p s» 0.35, we have p c ~ 1, and the 
critical disorder decreases as the particle density grows. It 
seems to approach zero as p becomes closer to 1/2. 

For fixed p, the average time T grows not only as 
the disorder increases towards p c , but also as the dis- 
order decreases and approaches zero. This implies that, 
remarkably, there is an intermediate value of p for which 
T is minimal. This value depends on the density, and de- 
creases as p grows. At the same time, the minimum be- 
comes sharper. For small densities, the curve is so flat 
that the minimum cannot be discerned, due to numeri- 
cal fluctuations in the value of T as \i varies. For larger 
densities, on the other hand, the minimum is very well 
defined. Thus, at least for sufficiently high densities, the 
segregation process is fastest at an intermediate value of 
the jump probability. Neither pure diffusion nor random 
long-range jumps is more efficient, to drive the system to 
the fully non-frustrated state, than this optimal combina- 
tion of the two processes. 

Before proceeding further with numerical simulations, 
we turn the attention to the analytical study of our sys- 
tem. Our analytical approach makes it possible to explain 
-or, at least, to clarify to a certain extent- the numerical 
result obtained so far. Moreover, it provides a clue to the 
steps to follow in further simulations. 
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P w-p P W + W~W + P P 



Fig. 5. Pair-antipair representation of a given configuration 
on a 10-site domain. 



3 Analytical results 

3.1 Pair-antipair representation 

We have already introduced particle pairs as an alternative 
representation of frustrated particles. In this section, we 
use an extension of this representation as a convenient tool 
for the analytical treatment of our system. 

A particle pair has been assigned to each occupied site 
whose nearest-neighbour site to the right is also occupied. 
The number of pairs P is thus given by the number of 
frustrated particles with a neighbour particle to the right. 
The density of pairs is p = P/ L. In the same way, we assign 
an antipair to each empty site whose nearest-neighbour 
site to the right is also empty. The number of antipairs is 
denoted by P, and their density is p — P/L. 

To complete the representation in terms of pair-like el- 
ements, we need to introduce the number W~ of occupied 
sites whose nearest-neighbours to the right are empty, and 
the number W + of empty sites with occupied sites to the 
right. However, due to the presence of periodic boundary 
conditions, we necessarily have W~ = W + . It is therefore 
sufficient to consider their sum, W = W~ + W + , which 
stands for the total number of "walls" separating pair and 
antipair domains. The corresponding density is w = W/L; 
we also introduce q = w/2 = W~ /L = W + /L. 

Figure |5] illustrates the pair-antipair representation of 
a given configuration of the system. Note that, due to 
the fact that two contiguous elements of this representa- 
tion share a common site, not all configurations are pos- 
sible. Specifically, only eight combinations of contiguous 
elements can occur: PP, PW~, W~W + , W~P, W+P, 
W + W~, PW+, and PP. They correspond to the eight 
different three-site neighbourhoods in the original repre- 
sentation. All of them are represented in Fig. \5\ 

Expressing the total number of particles and of empty 
sites in terms of the numbers P, P, and W, it is possible 
to verify that the densities p, p, and q, and the particle 
density p are related according to 



p_ + q = p, 

p + q = 1 -p. 
These relations make it possible to write 

q = p-p, 

p = 1 + p - 2p. 



(1) 



(2) 



Consequently, the only independent quantities in the pair- 
antipair representation are the particle density p -which, 
we recall, is one of our parameters- and the pair density 
p. In the frame of this representation, we analyze now the 
special cases // = and p = 1. 
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P \P W~P\ W + P W~P P 
P.^W~W*W~W + P W~P P 
W~W*P W~W*P W~P P 
W~W*P P W~W + W~P P 
P W~W + P.W~W + W~P P 

p w~w + w~mp w~F\ P 



P^W~W + W~W + W~W + W~P 
WW*P^ W~W + W~W + W~P 



rvWrff f w~p\ 
w-w + w~w + w~w + w~w + w 



Fig. 6. A realization of the segregation process of five particles, 
with p = (pure diffusion), within a 10-site spatial interval. 
Time runs downwards. At each step, the moving particle is 
shown in gray. To the right, the corresponding pair-antipair 
representation is given. Boxes indicate pair-antipair annihila- 
tion events. Lines stand for pair trajectories. 



3.2 Pure diffusion: p = 

For p = 0, frustrated particles can only move to nearest- 
neighbour sites, in the case they are empty, and long-range 
jumps are not possible. The motion is therefore purely 
diffusive, with the constraint imposed by the exclusion 
principle. 

Figure illustrates a particular realization, during 11 
time steps, in a 10-site portion of a larger system. Five 
particles, all of which are initially frustrated, reach a fully 
non-frustrated state. To the left of the plot, we show the 
corresponding (9-site) pair-antipair representation. Pair 
trajectories are shown with lines, and pair-antipair an- 
nihilation events arc indicated by boxes. 

An exhaustive analysis of the possible configurations 
in small neighbourhoods reveals that the are three ele- 
mentary events able to change the configuration. In the 
pair-antipair representation, they are 



PW~W+ 

pw-p -> 

PW+P -> 



W~W+P, 

w-w+w-, 

W+W-W+. 



The first one can occur in both directions, and stands for 
the diffusive motion of a pair P. Note that P moves two 
sites to the right or to the left at each step. Antipairs P, on 
the other hand, are immobile. The second and third events 
stand for pair-antipair annihilation. They are equivalent, 
in the sense that they can be transformed into each other 
by spatial inversion. 

Thus, in the limit of pure diffusion our population is 
analogous to a reaction-diffusion system of two species, P 
and P. While P diffuses through short-range jumps, P re- 
mains immobile. Moreover, if P and P become sufficiently 
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close to each other, they may undergo binary annihilation: 
P + P^O. 

Note that the possibility of this reaction was already im- 
plicit in the second of Eqs. (J5J), which states that the dif- 
ference p — p is a constant of motion (determined by the 
particle density p). 

The analogy with a two-species reaction-diffusion sys- 
tem makes it possible to apply well-know results of dif- 
fusion-controlled bimolecular reactions |2] . Specifically, we 
know that if the densities of pairs and antipairs are iden- 
tical, p — p, the kinetics is anomalous. According to the 
second of Eqs. J2J), such condition holds for p = 1/2, where 
the only stationary state for our system corresponds to a 
periodic array of alternating occupied and empty sites. 
Under this condition, fluctuations dominate the kinetics 
0. Pairs and antipairs segregate into separated spatial 
domains and, from then on, annihilation takes place at 
the domain boundaries only. Thus, the kinetics is drasti- 
cally slowed down and the pair density decreases, for long 
times, as 

p~r^\ (3) 

Figure |3 shows that numerical results are in good agree- 
ment with this prediction. 

For p p the approach to stationarity is exponential. 
If p < p (p < 1/2) all pairs are annihilated, while a num- 
ber of immobile antipairs subsists. Hence, the population 
is able to reach a fully non-frustrated state. If, on the 
other hand p > p (p > 1/2), moving pairs survive while 
all antipairs disappear. The remaining density of pairs is 
p* = 2p — 1, because p — 0. Asymptotically, the popula- 
tion is found in a stationary but dynamical state with a 
non- vanishing fraction of frustrated particles, as shown in 
Fig.0 

The analogy with binary annihilation makes it also 
possible to predict, at least in qualitative terms, the time 
behaviour of the pair density for p > 0, as compared 
with the limit p, = 0. As long as the disorder remains 
close to zero, pure diffusion will dominate the transport 
of particles and, thus, pair and antipair domains will form 
as a consequence of the mutual annihilation of these two 
species. While the ensuing annihilation events will mainly 
occur at the domain boundaries, rare long-range jumps 
of frustrated particles will give origin to non-local events 
where distant pairs and antipairs may disappear simulta- 
neously (see also Sect. 13. 3|) . These non-local events take 
place when a frustrated particle, jumping from a domain 
occupied by pairs, lands at a antipair domain. Their con- 
tribution to the variation of the pair density p is pro- 
portional to the product pp of the densities of pairs and 
antipairs, and to the disorder p. For p = 1/2 (p = p), 
this contribution implies an algebraic decay p ~ f ; for 
p / 1/2 the decay is exponential. In both cases, it accel- 
erates the decay of the pair density and, thus, the asymp- 
totic state is reached faster. Therefore, the evolution time 
T must be smaller for p > than for p — 0, as it was 
shown to be the case in our numerical simulations (Fig. 
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Fig. 7. All the long-range jump events which create or destroy 
particle pairs (up to spatial-inversion transformations). The 
moving particle is shown in gray. 



3.3 Full disorder: p = 1 

If p = 1, frustrated particles jump to randomly chosen 
empty sites, and the small- world disorder is maximal. At 
each time step, a particle can move a distance of the order 
of L. Regarding the segregation process, however, these 
random jumps represent an inefficient way of reaching a 
fully non-frustrated state. As our preliminary numerical 
simulations have shown, this results into a substantial in- 
crease of the times needed to accomplish full segregation. 

Since jumps take place to any empty site, it is possible 
to write an equation for the evolution of the pair density 
p by evaluating the probability of the events which create 
or destroy pairs in the limit of an infinitely large system. 
This approximation requires assuming that the positions 
of pairs and antipairs are not correlated. 

Figure shows -up to the transformation of spatial 
inversion- all the possible events with net creation or de- 
struction of pairs. The first line displays the only situa- 
tion where a pair -and the corresponding antipair- can 
be created. The moving particle starts from a site with 
one occupied nearest- neighbour (PW~) and lands at the 
empty site between two occupied sites (W~W + ). The fi- 
nal configuration at the original site is W~P, and two 
pairs (PP) have appeared at the landing neighbourhood. 
We represent this transition as 



pw-\w~w + 



W'PIPP, 



where to the left (right) of the arrow we indicate the initial 
(final) configurations of both the starting and the landing 
neighbourhoods. The (positive) contribution of this tran- 
sition to the average variation of the number of pairs P 
during a time step At is proportional to the probability 
of finding the initial neighbourhood configurations PW~ 
and W~W + when choosing a particle at random. As for 
PW~ , the probability of finding these two elements in con- 
tiguous positions -i.e., sharing a particle- is pq/p 2 - Equiv- 
alently, the probability of finding W~W + when choosing 
an empty site for landing, is q 2 /(l — p) 2 ■ Thus, the con- 
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tribution of the above transition to the variation of P is 

,3 



pq 



(4) 



P 2 (l-P) 2 

Spatial inversion transforms this transition into 

W + P\W~W + -> PW + \PP. 

By symmetry, its contribution to the variation of the num- 
ber of pairs is the same as above, so that APi will enter 
the evolution equation for the pair density multiplied by 
two. 

The three last lines in Fig. [7| stand for events where 
net destruction of pairs -and antipairs- occurs. They cor- 
respond to the transitions 

pw~\pp -> w-p\w+w-, 
pp\w-p -> w-w+\pw-, 
PP\PP -> w~w + \w + w-, 

respectively. In the first two of these transitions only one 
pair is destroyed, while in the third transition two pairs are 
simultaneously destroyed. Following the same arguments 
as above, we can calculate the respective (negative) con- 
tributions to the time derivative of the pair density, as 



AP 2 = -pp 2 q/p 2 (l-p) 2 , 
AP 3 = -p 2 pq/p 2 (l-p) 2 , 
AP 4 = -2p 2 p 2 /p 2 (l -pf 



(5) 



The contributions AP 2 and AP3 must be multiplied by 
two, to take into account the transitions obtained by spa- 
tial inversion. On the other hand, the transition giving 
origin to AP 4 is invariant under such transformation and, 
therefore, there is no further contribution associated with 
it. 

In summary, within our approximation, the average 
variation of the number of pairs P during a time step is 
given by AP = 2(APi + AP 2 + AP 3 ) + AP 4 . The time 
derivative of the pair density p = P/L can be evaluated as 
p = AP/LAt = pAP, where we have taken into account 
that At = l/N. Using Eqs. ©, ©, and ©, we get 



2p 



P(l-P) 



— — (3p 2 -3pp- p+p 2 ). 



(6) 



This equation has three fixed points, corresponding to sta- 
tionary values of the pair density. They are 



P± 



P ± \ P~ lP 



(7) 



Linear stability analysis shows that p* + is unstable for any 
value of the pair density p. For < p < 1/3, p*_ is neg- 
ative and unstable, while p$ is stable. For p > 1/3, these 
two fixed points interchange stability with each other, and 



p*_ > 0. Therefore, we have a transcritical bifurcation at 
p c = 1/3 ra 0.33, in good agreement with our numeri- 
cal results. If the particle density is lower than p c , the 
pair density vanishes asymptotically, and the system is 
able to reach a state where no frustrated particles subsist. 
For larger particle densities, on the other hand, a fraction 
p*_ > of the population remains always frustrated. It is 
remarkable that, according to our prediction, an infinitely 
large population is unable to find a fully non-frustrated 
state for 1/3 < p < 1/2, even when such states do ex- 
ist within that density interval. This feature makes clear 
that too much disorder in transport is inefficient as for the 
achievement of full segregation. 



3.4 Summary of analytical results 

It is useful to briefly summarize our analytical conclusions, 
in order to clarify to which extent our previous numerical 
results have been explained and to define the aims of fur- 
ther simulations, (i) For p = 0, we have established the 
equivalence between our system and the mutual annihila- 
tion of two diffusing species. Known results for the latter 
make it possible to explain the observed algebraic decay of 
the pair density, p ~ i -1 / 4 , for p = 1/2 (Fig. We have 
also shown that, for p > 1/2, the asymptotic pair density 
should be p* = 2p — 1. (ii) For p > 0, it was possible to 
show qualitatively that the evolution to full segregation is 
faster than for p = 0, as already found in the simulations 
(Fig. |2J|. (iii) For p = 1, we have derived an evolution 
equation for the pair density p. This equation makes it 
possible to predict, in particular, the asymptotic values of 
p as a function of p, from Eq. Q. Also, it predicts a tran- 
scritical bifurcation between the regimes of low and high 
particle densities at p c = 1/3. Combining conclusions (ii) 
and (iii), we can give analytical support to the observation 
of a minimum in the average time T as a function of the 
disorder p (Fig.HJ. In fact, while T decreases for small p, 
it diverges at a finite disorder level, p c < 1. 



4 Further numerical results and phase 
diagram 

Following our analytical results, it is first pertinent to 
study in more detail the asymptotic pair density p* for 
different values of the small- world disorder p and of the 
particle density p. Figure[S]shows numerical measurements 
of p* as a function of p, for four values of p. For p = 0, 
results are in excellent agreement with the analytical pre- 
diction, p* — 2/5—1, shown as a line. As p grows, the 
dependence of p* on p becomes increasingly non-linear 
and, at the same time, the critical density p c shifts to the 
left. For p — 1, the numerical value of the critical den- 
sity is slightly above the prediction p c = 1/3. Moreover, 
while analytical and numerical results for p* as a function 
of p are in good agreement, the analytical prediction is 
systematically larger. As shown in the inset of the fourth 
panel of Fig. [SI this discrepancy is more important close 
to the critical point. 
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Fig. 8. The stationary pair density p* as a function of the 
particle density p, for four values of the small-wolrd disorder 
p. For p = and 1, full lines show the analytical predictions 
(see Sects. |3~21 and HOI . For p = 1, the inset shows a close-up 
of the main plot for 0.3 < p < 0.5. 
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Fig. 9. The evolution of the pair density p for various values 
of the particle density p, with p — 1. Full lines correspond 
to simulation results, and dotted lines stand for the numerical 
solution of Eq. JBJ. 



As may be expected, the same effect is present in the 
comparison between the evolution of the pair density for 
p = 1 as obtained from simulations, and the solution of the 
evolution equation ©. Figure [5] shows, as full lines, sim- 
ulation results for various values of the particle density p, 
in the same conditions as in Fig. EI Dotted lines stand the 
corresponding solutions of Eq. ©. The latter have been 
obtained through a fourth-order Runge-Kutta integration 
of the equation with the initial condition p(0) = p 2 , which 
corresponds to a random distribution of particles. The 
agreement between the two results is very good for small 
densities, well below the critical value p c — 1/3. As p ap- 
proaches p c , however, a noticeable discrepancy arises in 
the long-time evolution of the pair density. On the other 
hand, the initial evolution obtained from simulations is 
well described by the analytical prediction. The long-time 
discrepancy seems to reach its largest relative value pre- 
cisely at the critical density p c . From then on, its effect 
decreases steadily, though a difference in the asymptotic 
values p* persist, as already shown in Fig. EI 

We identify two possible causes for the discrepancy be- 
tween simulation and analytical results for p, = 1. First, 
simulations are necessarily affected by finite-size effects. 
In our system, they convey the possibility that, due to 
fluctuations, the population reaches a fully non- frustrated 
state faster than predicted by the average-like arguments 
of the analytical formulation, which hold for an infinitely 
large system. We have already stated that, indeed, a fi- 
nite system should always be able to reach a fully non- 
frustrated state in a finite time for any particle density 
lower than 1/2. Such "early" equilibria may explain the 
fact that simulation results for the pair density are sys- 
tematically below the analytical prediction. As expected, 
the effect is relatively larger close to the critical point. 

The second source of discrepancy lies on the uncon- 
trolled assumption made to formulate Eq. © that the po- 



sitions of pair-like objects are not correlated. While long- 
range correlations are surely negligible for the random ini- 
tial condition, the evolution -which systematically tends 
to destroy frustrated configurations and to stabilize non- 
frustrated states- is expected to create such correlations. 
This is illustrated in Fig. El (for p — 0), where the creation 
of long-range structures such as • • • W + W~W + W~ • • • is 
apparent. While the effect of reaching "early" equilibria, 
discussed above, disappears for p > 1/2, the effect of cre- 
ation of correlations between pair-like objects is present 
in the whole domain of particle densities. 

Our main results for segregation in annealed small 
worlds are summarized in the phase diagram of Fig. 1101 
The parameter space (p,p), spanned by the small- world 
disorder and the particle density, is divided into two main 
zones. For low densities, we have the phase where the pop- 
ulation is able to reach a fully non-frustrated state while, 
for large densities, some degree of frustration persists even 
at asymptotically long times. Within the former, we find 
the curve over which the average time needed to reach 
the fully non-frustrated state is minimal. For a given par- 
ticle density, it gives the optimal small-world disorder to 
achieve segregation. In the frustrated phase, Fig.llOlshows 
curves corresponding to constant values for the density of 
pairs persisting at long times. 



5 Conclusion 

In this article, we have presented and studied a model 
for microscopic segregation in a homogeneous population 
of particles moving on a one-dimensional annealed small- 
world network. Motion is thus a mixture of short-range 
diffusion and long-range jumps, picking (frustrated) par- 
ticles with occupied neighbours and transporting them to 
empty sites. Small-world disorder is here identified with 
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Fig. 10. Phase diagram in parameter space. The bold line sep- 
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from that where frustrated particles subsist for asymptotically 
long times. The dotted line indicates the values of (fi, p) for 
which the fully non- frustrated state is reached fastest. Lines 
in the frustrated phase are level curves of constant stationary 
pair density p*. 



the probability of long-range jumps. While no macroscopic 
separation of phases takes place, for low densities the pop- 
ulation is driven toward a (fully non-frustrated) spatial 
distribution where particles are mutually separated by at 
least one empty site. 

From numerical simulations of the system we have ob- 
tained, on the one hand, the time needed to achieve full 
segregation in the regime of low densities. On the other 
hand, for high densities, we have analyzed the asymptotic 
number of persisting frustrated particles. In the parameter 
space spanned by the small-world disorder and the parti- 
cle density, we have determined the boundary between 
the two regimes, which interrelates the two parameters 
in a non-trivial way. These numerical results have been 
partially explained by means of an analytical approxima- 
tion, based on a representation of the particle distribution 
through pair-like elements. In this representation, the mo- 
tion of individual particles is associated with events of pair 
motion or pair-antipair annihilation. The approximation 
is useful for vanishing small-world disorder, where the seg- 
regation process becomes equivalent to the binary annihi- 
lation of two diffusing species. In this limit, our approx- 
imation makes it possible to predict the anomalous time 
decay of the number of frustrated pairs at the bound- 
ary between the regimes of low and high densities. In 
the opposite limit of maximal small-world disorder, where 
particle motion occurs purely through long-range jumps, 
the approximation predicts a transcritical bifurcation at 
the boundary. Numerical measurements of the asymptotic 
density of frustrated pairs are in good agreement with the 
analytical prediction. Also, the short-time evolution of the 



pair density is reasonably well approximated by the ana- 
lytical results. 

Our most remarkable results concern the regime of 
low densities. As stated above, in this regime, the aver- 
age time needed to reach a fully non-frustrated state is 
a meaningful overall quantity characterizing the dynam- 
ics. We have found, first, that this time attains a mini- 
mum at an intermediate value of the small- world disorder. 
In other words, there is an optimal combination of dif- 
fusion and long-range jumps that makes the segregation 
process fastest. Second, both numerical and analytical ev- 
idence show that, in the limit of an infinitely large system, 
there is a critical value of the small-world disorder above 
which the time needed to reach the fully non-frustrated 
state diverges. Too many long-range jumps, thus, make 
the achievement of full segregation impossible. The critical 
value of the small- world disorder determines the boundary 
between the regimes of low and high densities. 

Besides the relevance of these conclusions as for the 
dynamics of segregation itself, they are particularly inter- 
esting from the viewpoint of small- world based processes. 
In fact, it is well known that the threshold of critical be- 
haviour in a large class of systems with underlying small- 
world structure shifts toward zero small-world disorder 
as the system size N grows [12] . This is a direct conse- 
quence of the fact that the cross-over from the ordered-like 
regime to the random-like regime in the geometrical prop- 
erties of small- world networks occurs at a disorder level of 
order TV -1 |13II14) . Only a handful of small-world based 
processes have been reported where critical behaviour or 
non-monotonic dependence on the disorder persist at fi- 
nite disorder levels as N is increased [51l51ll5 | . Our system 
is a novel instance of this rare category, combining both 
non-monotonic dependence and critical behaviour. 

Several generalizations of the present model would be 
worth considering. First, the same dynamical rules can 
be implemented in two- or higher-dimensional lattices. In 
more dimensions, an additional degree of freedom is given 
by the various ways in which the neighbourhood of a site 
can be defined. In turn, this determines several possibili- 
ties as for the definition of frustrated and non-frustrated 
states. Another important generalization has to do with 
the inclusion of noise. In the present version of the model, 
diffusion events are permitted only when they drive a par- 
ticle from a frustrated state to a non-frustrated state. In 
this sense, diffusion in our model is unidirectional toward 
segregation. Noise can be implemented by allowing a cer- 
tain degree of reversibility for diffusion events, such that 
non-frustrated particles might occasionally be lead to frus- 
tration. This additional mechanism would be particularly 
relevant for small disorder levels, when motion is dom- 
inated by diffusion. In the pair-antipair representation, 
reversible diffusion events would be represented by spon- 
taneous creation of a pair and an antipair, — > P + P, 
corresponding to the transition W~W + W~ — * PW~P 
or to its spatially inverted form W + W~W + — ► PW + P. 
Finally, we recall from the introduction that segregation 
can be thought of as a process aimed at distributing a 
population in order to optimize the exploitation of space. 
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Our model may be extended to drive the population to a 
state where the inter-particle distance is maximal, under 
the constraint imposed by the particle density. 
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